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COMPRESSIBLE-FLOW' SOLUTIONS FOR THE ACTUATOR DISK 
By James B. Delano and John L. Crigler 


SUMMARY 


Generalized solutions for the actuator disk in subsonic compress- 
ible flow have been obtained and the compressible- and incompressible- 
flow phenomena are compared. The results show large differences in the 
flow variables across the disk. 

Two primary types of solutions were obtained: one wherein the flow 

Imme diately behind the disk is subsonic and the other wherein this flow 
is supersonic. For the subsonic wake, the maximum efficiencies for com- 
pressible and incompressible flow are the sams up to choked flow. The 
maximum efficiencies for the supersonic wake are lower than for the sub- 
sonic wake but this difference decreases as the free -stream Mach n umb er 
approaches 1.0. 

The use of actuator disks in tandem is one method of delaying the 
large losses in efficiency associated with choked flow, especially at 
low values of Mach number and high loadings. A level of efficiency 
almost as high as for incompressible flow could be maintaine d until the 
flow in all disks became choked. A single actuator disk gives the ideal 
efficiency only for power loadings up to those required to choke the flow. 


INTRODUCTION 


The simple momentum theory has been used to great advantage in 
analyzing propeller problems and predicting propeller efficiencies for 
known operating conditions. In the simple momentum theory the propeller 
hag been treated as an actuator disk in incompressible flow with the 
thrust of the propeller uniformly distributed over the disk and with a 
constant axial interference velocity over the disk. Since, in practice, 
propellers must operate in compressible flow, it was felt that an exten- 
sion of the actuator -disk theory to compressible flow might prove useful. 
Consequently, the p imp le momentum theory for an actuator disk has been 
extended to compressible flow based on what appear to be the logical 
assumptions . 
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If a large pressure change is assumed to occur across the disk and 
this change occurs isentropically, then a large density change must also 
occur. Under these conditions, obviously, an increase in density across 
the disk must also be accompanied by a decrease in velocity through the 
disk to satisfy the conservation of mass flow. Solutions for the actu- 
ator disk in compressible flow are presented in reference 1 for one Mach 
number and for a limited range of power loading. He study of refer- 
ence 1 did not include the possibility that the flow anywhere in the 
stream have a Mach number greater than 1. 

A more comprehensive study of the flow phenomena and performance 
characteristics for an actuator disk in compressible flow is presented 
herein. He scope of the treatment has been extended to include the 
effect of entropy changes associated with profile and shock losses. In 
addition, a study is made of the conditions under which the flow in the 
slipstream exceeds a Mach number of 1.0, and the performance of the disk 
is calculated. 
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SYMBOIS 

(y p\l/2 

speed of sound, (— J > ft/sec 

slipstream cross-sectional area 
total energy, ft-lb/slug 
acceleration of gravity 
total pressure, Ib/sq ft 
mass flow, slug/sec 
Mach number in slipstream 
Mach number behind normal shock 

static pressure, lb/sq ft 

static -pres sure rise through disk, p^ - p^, lb/sq ft 

total power added to slipstream by actuator disk, ft-lb/sec 
power producing thrust, ft-lb/sec 
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power expended in profile or shock losses in the disk, 
ft-lb/sec 


power disk-loading coefficient. 


power disk-loading coefficient. 


power disk-loading coefficient, 


Wo\ 

p i 

K t o 3a x 


value of P c when the flow becomes choked 


thrust, lb 

thrust disk-loading coefficient. 


velocity in slipstream, ft/sec 

ratio of specific heats (for air, 1. It-) 

air density, slugs/cu ft 

TV 0 T c 

efficiency, = — 

p p c 

ideal induced efficiency for compressible flow 
ideal induced efficiency for incompressible flow 

H 1 - 


from P Ci = 


^i 


inc 


Subscripts : 

0 station far ahead of disk 

1 station immediately ahead of disk 


k 
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2 station immediately behind disk 

3 station far behind disk 

F conditions pertaining to front disk of disks in tandem 

R conditions pertaining to rear- disk of disks in tandem 

s solution for supersonic wake 

THEORY 

Basic Assumptions and Considerations 


The actuator disk and the stations used in the analysis are shown 
in figure 1. The assumptions made in reference 1 are used herein. These 
assumptions are: 

(a) The propeller may be represented by an actuator disk of zero 
thickness. 

(b ) The velocity distribution over any cross section of the stream- 
tube passing through the disk is taken to be uniform (one -dimensional 
flow ) . 

(c ) The energy is added to the slipstream instantaneously and 
uniformly over the disk. 

(d) The slipstream area is continuous through the disk so that the 
cross-sectional areas Immediately ahead of and immediately behind the 
disk are equal. 

(e ) The actuator disk introduces no rotation in the slipstream. 

(f ) In the final wake the static pressure returns to the free- 
stream value = Pq)* 

(g) The flow ahead of the disk is isentropic. 

The scope of the treatment has been extended to include solutions 
wherein arbitrary energy losses , such as might be associated with blade 
profile drag, occur at the actuator disk. The solutions presented herein 
are for free -stream subsonic flows wherein the Mach number immediately 
ahead of the disk does not exceed 1.0. In order for M-j_ to reach values 

greater than 1.0, the slipstream must contract and then expand ahead of 
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the disk (similar to the shape of a Laval nozzle) with a Mach number 
of 1.0 at the minimum cross section. This type of flow seemed hi ghly 
unlikely for a free boundary even though pressure waves from the disk 
could possibly travel upstream outside the slipstream boundary. 

The addition of power greater than that required to produce 
Ml = 1.0 is assumed, therefore, to cause no further changes in the 

flow upstream of the disk. With this assumption, solutions of the 
problem consistent with the preceding assumptions are possible only 
for nonisentropic flow between stations 1 and 3 . Under such condi- 
tions, depending upon the assumed distribution of the losses between 
stations 1 and 3 there are two possibilities: (a) the Mach number Mg 

Immediately behind the disk (station 2 , see. fig. l) may decrease or 
(b) it may become supersonic. Hie possibility of a supersonic wake has 
received some consideration in reference 2 . 

Solutions have been obtained under the assumption of both subsonic 
and supersonic exit velocities, and the power disk loading of the prob- 
able change from subsonic flow to supersonic flow behind the disk, for 
a given stream Mach number, has been determined. 


Method of Solution and Basic Equations 

Subsonic wake . - Figure 1 defines the four stations in the slip- 
stream studied in the course of -the analysis. Free -stream conditions, 
which are known, exist at station 0. Stations 1 and 2, respectively, 
are immediately ahead of and immediately behind the disk. The power is 
added to the slipstream between stations 1 and 2. Station 3 is in the 
final wake where the static pressure is again equal to the free -stream 
static pressure. The problem is to determine the flow variables through- 
out the slipstream associated with given free-stream conditions and power 
disk loading. 

Hie solution for a given set of operating conditions is obtained 
when the flow variables satisfy the thrust, energy, and continuity 
relationships. Since the relationships of the various flow variables 
are not explicit, a method of successive approximations is employed to 
effect a solution. Hie following method of solution was employed. 

A tentative assumption of Mach number at station 1 permits the 
determination of all the flow variables at station 1 from isentropic 
considerations. Hie variation of streamtube cross-sectional area with 
Mach number for isentropic flow is given in reference 3 as 
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Aq " KlI 


1 + 


7-1 


7+1 

2(7-1) 


7-1 0 

1 + — - — Mq 2 


(1) 


The value of is never assumed to he greater than 1.0 for the final 

results presented in this paper for the reasons previously stated. The 
static pressure p^ is obtained from the isentropic relations between 

the total and static pressures between station 0 and 1 given by 


and. 



(2) 


Ho - % 


(3) 


The density at station 1 is obtained from the isentropic pressure- 
density relation 



00 


The velocity is obtained from the continuity of mass flow 


h. PqAd 
v o " p x A i 


(5) 


The total energy at station 0 is equal to the total energy at station 1 
and is given by 

E 0 =E 1 1 
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Hie flow variables at station 3 are found next. These variables 
depend on the amount of energy being added to the slipstream at the 
disk, how this energy is added (isentropic or nonisentropic compression), 
and the mass flow. The total energy at station 3 is equal to that at 
station 1 plus the energy added by the disk. 

P 

Since is the power per unit mass flow added to the stream 

"A’i 

"by the disk.. 


V 0 2 + 7 fo + F _ 7 3 2 , 7 ?3 

2 7 - 1 P 0 Pi'S. 7 ! 2 7 ~ 1 P 3 

This equation can be rearranged to give 

P = v 3 2 ~ Y 0 2 + y p o/ P Q _ A 

P 1 A 1 V 1 2 y-l P^\P 3 “ J 

The total power can be divided into two parts 

P = Pl + P 0 


(7) 


( 8 ) 


where 


pi A i v i 



(9) 


is the power per unit mass flow required to produce thrust and 

p Q = _r_ £o /fo . A 

p l A l v l 7 " 1 p 0 \ p 3 / 


( 10 ) 


is the additional power associated with entropy changes resulting from 
profile or shock losses in the disk. 

For the case where no entropy change takes place, = Pq and 
hence P Q = 0. Then, only induced losses are incurred and P = P.j_. The 
velocity at station 3 can be obtained from equation (9). When Pi is 
changed to the power disk-loading coefficient P c _^, the velocity ratio 


becomes 
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v o 



(H) 


Sijnilarly, the density ratio — 

p 3 

which can he put in the following 


can he obtained from equation (10) 
form: 


^0 

p 3 


1 + 


(y - i) 

2 



F °o 

A<)/ A l 


( 12 ) 


also 


P 


c 



The speed of sound at station 3 is given by 


a 3 “ 



and the Mach number is given by 


Mo 



(13) 


(lA) 


(15) 


The total pressure is found from 



(16) 


where P 3 = p Q . 


The area ratio is given "by 


f3 = P 0 V 0 

Ao “ P3V3 


(17) 


All' the flow variables at station 3 have now been determined for 
the assumed Mjj_ and for the given power per unit mass flow, that is, 

P c is constant. 
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Conditions at station 2 are next determined by isentropic-f low- 
relations in working backwards from, station 3 to station 2. The area 
ratio A 3 /A 2 is given by 


A 3 _ A 3 Pl V l 
*2 " ~ 


(18) 


When isentropic flow occurs across the disk, equals Pq. The Mach 

number at station 2 is determined from the area ratio 


The static pressure 


h. 

P 2 


l 

M 3\i + M2 2 

is obtained from 


7+1 

2(7-1) 


7 



where 

% - h 3 

The density at station 2 is obtained from 



The velocity V 2 is obtained from the continuity of mass flow 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Vo = 


p l V l 


(23) 


since A-[_ = A 2 . The energy at station 2 is equal to the energy at sta- 
tion 3 and is given by 
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E, 


= Eo 


V- 


— + — 2 — _£ 

2 7 - 1 P 2 ^ 


_3 , y ?3 

2 7 - 1 p 3 


(24) 


The determination of the flow variables at station 2 is now complete 
for the assumed value of 


The corresponding value of thrust consistent with the assumed 
value of M^, P c , and P c ^ is found from the sum of the pressure 

forces and momentum changes across the disk 

T - »A V l( T 2 - V l) + - p l) A l < 2S > 

It can "be shown., however, “by methods similar to those shown In refer- 
ence 4 that the net pressure force on the fluid between stations 0 and 3 
must be zero. For this condition, the thrust is given also by the change 
in momentum between stations 0 and 3 

T - - V 0 ) (26) 

In general, the arbitrarily assumed value of M-j_ associated with a 

given value of P„ and P_ will not satisfy equation ( 2 6 ). The 
c ^o 

correct final solution will be obtained when Mj_ has been chosen in 

such a way that equations (25) and (26) give the same value of thrust. 

If the thrusts given hy these two independent thrust equations (25) 
and (26) are plotted against for given values of P c and P Cq , 

the intersections of the two curves will give the solution. 

The thrust loading T/A-^ obtained from the two independent thrust 
equations (25) and ( 26 ) is shown plotted in figure 2 against M-j. for 
constant values of F c , P Q /P, and free -stream. Mach number. Three 

typical cases are presented to show the types of solutions obtained. 

In figure 2(a) there are two values of (one subsonic and one super- 

sonic ) that satisfy the two independent thrust equations; hence, there 
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are two possible solutions. In line with, the previous discussion, only 
the subsonic values of M]_ will be used. In figure 2(b) there are no 

crossings of the two thrust curves and thus no solutions. Ihis fact 
indicates that there is a limit to the maximum power disk loading for 
which solutions can be obtained for isentropic flow throughout the slip- 
stream. For example, if a sufficient part of the power had gone into 
losses, a solution would have been possible, such as is shown in fig- 
ure 2(c). In figure 2(c) all conditions are the same as in figure 2(b) 
except that the power loss P Q /P (drag or shock loss) is 26.4 percent 

of the total power. In this case only one solution was obtained and 
that is for Mj_ * 1.0. If P Q /P is increased for the same values 

of P c and Mq, two solutions will be obtained, one with subsonic 

and one with Mq_ supersonic. 

These results can be summarized as follows: When P c is less 

than P c *, there are two values of M?j_ (one subsonic and one super- 
sonic ) that will satisfy the two independent thrust equations. When P c 
is equal to P c *, there is only one value of Mj_ that produces a solu- 
tion, that is, = 1.0, and the flow is everywhere isentropic. When 
P c is greater than P c *, no solutions are possible unless arbitrary 
losses are introduced in the disk. These losses are similar to losses 
associated with the blade profile drag and shocks. Consequently, the 
flow through the disk is no longer isentropic. 

Thus far only solutions having subsonic wakes have been considered. 
It will be shown later under "Results and Discussion" that there is a 
transition from a subsonic wake to a supersonic wake at some value of 
P c greater than P c *, that is, after the flow into the disk becomes 

choked. 

Supersonic wake .- In order to analyze the flow behind the disk 
with a supersonic wake, the flow out of the disk is assumed to be 
similar to the supersonic flow out of a nozzle in which the discharge 
or exit pressure is less than the ambient pressure at the exit. The 
attendant shock formations in the wake are discussed in detail in refer- 
ence 3 (pp. 172-174). This compression gives rise to shock losses in 
the wake so that the treatment used for the subsonic-wake solutions no 
longer is applicable. The following method has been used to obtain the 
supersonic -wake solutions. 
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All the flow factors through station 1 have been determined since 
the flow is choked. A value of V 2 (Mg > l.o) is assumed. The density 

at station 2 is obtained from the continuity of mass flow 


p 2 “ 


Pl V l 

v 2 


p 

In general, by assuming — , p„ can be obtained from the energy 

m u 

equation 

E 0 + 1 = 
or 



( 27 ) 


The speed of sound is given by 


®2 



and the Mach number by 

Vo 

“2 = 4 

The thrust is given by the sum of the pressure forces and momentum 
changes across the disk 


T 

m 


Vo - V, + 


P 2 ~_Jl 
p l V l 


which is another form of equation (25). The efficiency is given by 


n = 


2 

m 


V 0 


P 

m 
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The power disk-loading coefficient is given by 

P Aq • 



2 


The solution for obtaining all the flow variables across the disk, 
the efficiency, and power disk-loading coefficient is now complete. It 
is not necessary to make a detailed analysis of a balance of energy and 
thrust by comparing the total energy and thrust at stations 2 and 3* 

P 

The arbitrary assumption of and — produces an infinite num- 

Ct lU 

her of solutions which in general do not represent the maximum-efficiency 
case. The maximum efficiency, however, is of particular interest and 
occurs when the flow through the disk is isentropic. For this case, pg 

is obtained from the isentropic pressure-density relation 


P 2 



The value of — corresponding to this isentropic process can now he 

obtained directly from equation (27). It is to be noted that in the 
subsonic -wake case, there is a definite limit to the value of the power 
(P c = P c *) the disk can absorb isentropically. For the supersonic-wake 

case, however, much larger amounts of power can he absorbed isentropically 
by the disk. 


RESULTS AND DISCUSSION 
Ideal Solutions 


Solutions of the flow variables throughout the slipstream are pre- 
sented in figures 3 "to 8 for forward Mach numbers of 0.60, 0-70, 0.80, 
0.85, 0.90, and 0.95. These results are presented as Mach number and 
as ratios of pressure, density, velocity (referred to free-stream condi- 
tions), and ratio of cross-sectional area of the slipstream (referred to 
the area of the actuator disk), all as functions of. the power disk- 
loading coefficient P c . The variation of efficiency with P c is pre- 
sented in the same figures as the pressure ratios. When the flow into 
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the disk is subsonic (m^ < l), the flow "behind the disk (station 2) must 

he subsonic, as will be shown later, and only the solutions for a sub- 
sonic wake were obtained. When the flow into the disk becomes choked 
( = l), solutions were obtained for both subsonic and supersonic flows 

at station 2. For the subsonic -wake case, however, solutions can be 
obtained only by introducing losses in the disk. All curves for solu- 
tions with subsonic flow at station 2 are marked with plain symbols 
while the curves obtained for supersonic flow at station 2 are marked 
with symbols having subscript s (figs. 5 to 8). Choked flow occurs 
at the lower value of P c at which there is a sharp break in the curves 

There is a range of power disk-loading coefficient for which both 
subsonic and supersonic wakes appear as alternative possibilities. 

Since there is nothing in the equations to indicate which type of flow 
might occur and which efficiency might be obtainable, the results were 
studied to determine when the change from a subsonic wake to a super- 
sonic wake could be expected. A supersonic wake is possible for all 
exit velocities for which a shock pattern exists at the exit or extends 
downstream from the exit. The lowest supersonic -wake Mach number which 
can exist is seen by analogy with nozzle flow to correspond to the con- 
dition of a normal shock immediately behind the actuator disk or at the 
exit of a supersonic nozzle. In this case, except for an infinitesimal 
distance behind the actuator disk, the wake flow is actually subsonic. 
Consequently, the subsonic -wake calculations should be applicable to 
this case and the losses which it is necessary to assume in order to 
obtain a solution correspond to the losses in a normal shock at the 
exit. 


A direct calculation of the Mach number immediately behind the 
disk for the case with a normal shock at the exit is difficult to make 
because the boundary pressure at the exit is not known at this stage of 
the calculation. However, since the subsonic- and supersonic-wake calcu 
lations should give the same results (corresponding efficiencies and 
powers are equal) for the normal shock at the exit, the intersection of 
the efficiency curves for the two types of solution (subsonic- and 
supersonic -wake solutions) should correspond to the condition of a 
normal shock at the exit. In figure 9, the Mach number behind a normal 
shock is plotted against the Mach number in front of a normal shock (see 
data from table 3 of ref. 3 )j this relation is shown by the solid line. 
The abscissas of the individual points (shown by symbols ) which are 
plotted in figure 9 are the values of Mg from the subsonic-wake solu- 
tions at the point of intersection of the efficiency curves for the 
subsonic-wake and supersonic -wake solutions. The ordinates of these 
points are the values of (Mg)^ taken from the supersonic -wake solutions 

at the corresponding conditions. It is seen that the relation between 
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the subsonic and supersonic values of Mg at the point of intersection 
of the efficiency curves is exactly that for a normal shock. 

The curves, figures 5 to 8, for all flow parameters are shown as 
solid lines for the subsonic-wake solutions up to the transition value 
of P c and as dashed lines at higher values of P c . The solutions for 

the supersonic wake are shown as solid lines, beginning at the value 
of P c for equal efficiencies from the two wake solutions and mi higher 

values of P e . At lower values of P c , the solutions for the supersonic 
wake are shown as dashed lines and are presented only to show the range 
of calculations. 

The curves in figures 3 to 8 for compressible flow show two Impor- 
tant discontinuities of the flow parameters. The first discontinuities 
occur at the value of P c at which M-j_ = 1.0 which is called the 

critical power disk- loading coefficient and is designated as P c *. 

Choking cannot occur for the incompressible-flow case; therefore, the 
mass flow can increase indefinitely as P c is increased, whereas for 

compressible flow the mass flow reaches a maximum for a given forward 
Mach number when the Mach number immediately »b<*ari of the disk reaches 
a value of 1.0. The second break in the curves occurs at the value of 
P c where the efficiency for the subsonic-wake solution is equal to the 
efficiency for the supersonic-wake solution. 

Curves for the supersonic -wake solutions are not included for the 
lower Mach numbers (figs. 3 a^^d 4). In these cases extremely high 
power disk loadings are required to choke the flow and supersonic calcu- 
lations for choked flow were not presented as it was felt the power disk 
loadings were "beyond the practical range for propeller operation. 

Pressure . - The pressure rise through the disk for compressible and 
incompressible flow at power disk loadings below choking is compared at 
several forward Mach numbers in figure 10 for maxi mum-efficiency condi- 
tions (no profile losses). For a given value of P c and Mg, the pres- 
sure rise through the disk is much greater for compressible flow than 
for incompressible flow. The figure shows that the pressure rise for 

compressible flow is — — — times that for incompressible flow at 

(l - Mg 2 ) 

light power loadings. The (a) parts of figures 5 to 8 also show that, 
for a given forward Mach number when the flow into the disk reaches a 
Mach number of 1.0, practical ly no additional pressure rise can he 
obtained even for large increases in power for the subsonic -wake solu- 
tions. The same figures also show a large decrease in pressure across 
the disk for the supersonic -wake solutions. 
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Density . - The changes in density are similar to the changes in 
pressure . 

Velocity and Mach number .- The velocity changes are similar to the 
Mach number changes throughout the slipstream. The velocity in front 
of the disk is much greater than the incompressible -flow value and the 
difference increases with P c up to choking. (See fig. 6(c).) The 

most significant result is the great drop in velocity through the disk 
for the subsonic-wake solution, and the great increase in velocity for 
the supersonic -wake solution. The effect of forward Mach number and 
power disk loading on the velocity immediately behind the disk is pre- 
sented in figure 11 for the subsonic-wake case. This figure shows that 
the velocity Immediately behind the disk is less than the free-stream 
velocity for high subsonic Mach numbers for all power disk loadings. 

At lower free-stream Mach numbers (e.g. , Mg « 0.65 to 0 . 69 ) the velocity 

behind the disk is greater than the free-stream velocity. 

Figure 12 shows a comparison of the average -velocity ratio through 
the disk for incompressible and compressible flows up to choking. The 

V 1 + v 2 

average velocity through the disk 5 is higher for compressible 

flow than for incompressible flow, and increases with P c . At low power 
loadings, however, there is not much difference. 

From the great differences in the flow factors for incompressible 
and compressible flow the vise of conventional incompressible-flow theory 
for propellers might be expected to give erroneous results in calcu- 
lating propeller performance characteristics. Comparison of experimental 
and calculated results, however, for propellers even at high subsonic 
Mach numbers has shown remarkably good agreement in over-all character- 
istics for light loadings* This good agreement may be due to the fact 
that the average of and V 2 for compressible and incompressible 

flow is not greatly different for light power loadings. For high power 
loadings, however, experimental propeller characteristics may differ 
greatly from calculations based on incompressible -flow theory. 

Efficiency . - Figure 13 shows the variation of efficiency with power 
disk-loading coefficient for all free-stream Mach numbers for both the 
subsonic- and supersonic-wake solutions. The changes in efficiency for 
a given free-stream Mach number were discussed earlier in connection 
with transition of the flow from a subsonic to a supersonic wake. The 
drop in efficiency from the subsonic- to supersonic -wake case (tran- 
sition drop) becomes less as the free-stream Mach number increases. 
Consequently, at the high subsonic Mach numbers the efficiency level 
is quite high for the supersonic wake case. 
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In figure 13 one maximum efficiency point obtained in experiments 
in the Langley l6-foot transonic tunnel is included for a free -stream 
Mach number of O .96 and P c = 0.054 for comparison with actuator -dish 

results. This efficiency falls between the lines representing the 
subsonic-wake solution with the flow choked and the supersonic-wake 
solution. This efficiency is about l6 percent lower than the ideal 
efficiency calculated for the supersonic wake and this loss represents 
approximately the blade profile-drag loss. The rest of the losses would 
appear even without blade profile-drag losses and would show up as shock 
losses in the wake, lhe efficiency and power disk-loading coefficient 
obtained for this propeller show that this propeller has a supersonic 
wake immediately behind ‘the propeller plane. 


Effect of Additional Losses 

Subsonic wake .- Important effects are produced when additional 
losses are included in the analysis. These effects are shown in fig- 
ure 14 for a forward Mach number of 0.85 and for values of P Qi /P of 0, 

0.1, 0.2, and 0.3 for the subsonic -wake case. For a given value of P c 

the deviation of the flow factors from the free-stream values decreases 
as the profile losses are increased. For a fixed value of P o yp, the 
limit of the subsonic-wake solution occurs when the Mach number into 

the disk reaches a value of 1.0. This limiting condition represents 
choked flow. The locus of similar limiting points for all values of 
P Q /P is designated as the "choke line.” The choke lines are shown as 

dashed lines in figure 14. Once the flow becomes choked and the power 
is increased, solutions can only be obtained by the addition of further 
arbitrary losses in the disk, that is, choking occurs at higher values 
of P c . Transition of the flow from a subsonic to a supersonic wake 
must occur somewhere along the choke lines. 

The effect of additional losses on efficiency is shown in figure 15 
for forward Mach numbers from 0.70 to 0.95 for the subsonic wake. The 
variation of efficiency for constant values of P Q /P is shown by the 
dashed lines. It can be shown that the efficiency is given by 

T1 P n P n 

— =1 — where tu A is the value of T) when — = 0 for the 

Tip =Q P *0=° P 

value of P c under consideration. A s P c is increased above the value 
required to choke the flow, the efficiency decreases rapidly. 

Supersonic wake .- Similarly, for the supersonic -wake case, the 
addition of profile-losses results in lower efficiencies than shown in 
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figure 13 . The value of P c at which choking occurs will increase with 
an increase of profile loss and the transition of flow from a subsonic 
to a supersonic wake will occur at some point along the choke line at a 
higher value of P c than the transition for the zero-profile-loss case. 

Previously, it was pointed out that the loss in efficiency associ- 
ated with transition of flow from a subsonic wake to a supersonic wake 
decreases as the free -stream Mach number is increased. Consequently, 
it is evident that if sufficient power is added to produce a supersonic 
wake, even at a free -stream Mach number of 0.8, a very large loss in 
efficiency results (from 96 percent down to 70 percent, see fig. 13 ) . 

The problem of particular importance, then is the determination, if 
possible, of methods to increase the efficiency by delaying the 
adverse effects associated with choking the flow. This can be done 
by the use of actuator disks in tandem. 


Actuator Disks in Tandem 

Solutions for actuator disks in tandem require knowledge of the 
flow field between the disks. In order to simplify this problem, it 
can be assumed that the disks are widely spaced and that only the mass 
flow passing through the first disk passes through the following disks. 
In addition, it is assumed that the final wake conditions of the leading 
disk are taken as the free-stream conditions of the following disk. 

Such a configuration resolves the problem into the solution of a series 
of single disks for which the solutions already obtained for a single 
disk can be used. 

When a single actuator disk is operating at a condition where it 
is absorbing just enough power to produce choking, the addition of more 
power will produce a rapid decrease in efficiency. If, however, this 
additional power were absorbed by a second coaxial disk located far 
behind the first disk, it can be seen that the second disk could operate 
at a high level of efficiency as long as the flow into it were not 
choked. Consequently, the over-all efficiency of the two disks would 
be higher than for a single disk absorbing the same' total power. 

Solutions have been obtained for two, three, and an infinite num- 
ber of disks In tandem using the methods presented in the appendix. No 
attempt was made to present all the flow factors since only the effi- 
ciency that can be obtained at high values of total power is of interest 
in the following analysis. 

Calculations were made on the assumption that the flow through the 
leading disk has just choked. Additional power is then added; only to 
the following disk until it becomes choked. The calculations were 
terminated for a finite number of disks when the flow into the last 
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disk became choked and, for an infinite number of disks, when the flow 
of the final wake reached a Mach number of 1.0. In this process the 
addition of power to any disk does not alter the flow conditions at any 
preceding disk. Each added disk becomes smaller than the preceding 
disk. This method of analysis was used for mathematical expediency and 
does not affect the results. For the case being considered (Mq = 0.8o), 

the diameter of the second disk was reduced about 3 percent and that of 
the third disk about 4 percent. The results of this analysis are pre- 
sented in figure 1 6 for a forward Mach number of 0.80, as over -all 
efficiency against power disk-loading coefficient based on total-power, 
free-stream conditions for the front disk and the area of the front 
disk. The efficiency for tandem disks remains almost as high as the 
incompressible efficiency and is considerably higher than it is for the 
single disk absorbing the same total power. The efficiency for the 
t and em disks is lower than the incompressible value of efficiency 
because the mass flow into the tandem disks remains fixed as the power 
is allowed to increase. Choked flow for a single disk occurs at a value 
of P c = 0.174 and the further addition of power results in a rapid 
loss in efficiency until the flow changes from a subsonic to a super- 
sonic wake. For tandem disks this rapid loss in efficiency is delayed 
until the flow in all the disks is choked. Choked flow occurs at 
P c = 0.263 for two disks, P c = O .319 for three disks, and P c = 0.584 

for an infinite number of disks. Thus, a considerable increase in high- 
efficiency power range results from the addition of one disk; further 
increases in power range occur for additional disks but this gain 
decreases as the number of disks is increased. 

Since it was possible to obtain a high level of efficiency by 
adding the power in more than one disk, it is obvious that the effi- 
ciency obtained for a single actuator disk no longer represents the 
optimum or ideal efficiency for all operating conditions. 

The use of closely spaced disks rather than widely spaced disks in 
tandem represents a practical propeller configuration. Experimental 
results show that final wake conditions are reached a short distance 
behind propellers and it is believed the results obtained with widely 
spaced disks are indicative of what may be expected for closely spaced 
propellers. Tests are needed to determine the effect of spacing on the 
performance of tandem propellers. 


CONCLUDING REMARKS 


Comparison of the incompressible- and compressible-flow phenomena 
show large differences in the flow variables across the actuator disk. 
For the case where only isentropic flow is considered for an actuator 
disk, the maximum efficiencies for compressible and incompressible flow 
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were found to "be the same for power disk loadings up to those first 
required to choke the flow. 

When the flow becomes choked, the further addition of power pro- 
duces a rapid loss in efficiency. This loss in efficiency is associated 
with the transition of flow immediately behind the disk from subsonic to 
supersonic flow. When the flow Immediately behind the disk is super- 
sonic, the rate of efficiency loss with power disk loading is similar 
to but greater than that for incompressible flow. 

The use of actuator disks in tandem is one method of delaying the 
large losses in efficiency associated with choked flow, especially at 
low values of Mach number and high loadings. A level of efficiency 
almost as high as for incompressible flow could be maintained until the 
flow in all disks became choked- It is obvious, therefore, that a 
single actuator disk gives the ideal efficiency only for power loadings 
up to those required to choke the flow. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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APPENDIX 

SOLUTION OF SPECIAL CASE OF ACTUATOR DISKS IN TANDEM 


For coaxial actuator disks in tandem any two successive disks are 
ass ume d to tie far eno ugh apart so that the final, wake conditions of the 
front disk may he assumed to he the free-stream. conditions of the fol- 
lowing disk. Hie problem is simplified hy the assumption that there are 
no profile or shock losses incurred. In addition, only the mass flow 
passing through the first disk is assumed to pass through all the disks. 
This assumption requires that the area of the second disk he smaller 
than the area of the first disk and just the right size to work on all 
the mass flow. These assumptions make possible the use of solutions 
already presented for the single actuator disk. Considerations similar 
to those of reference 4 show that the form of the general momentum equa- 
tions Involving thrust and boundary forces for this case remains the 
same and that the net boundary force is zero when the static pressures 
far in front of and far behind the acturator disks are the same. Fig- 
ure 17 is a sketch of the stations used In this analysis. 


The over -all thrust of the system is given by 

T * m ( v 3R - v °f) 


(1) 


where 


is the final wake velocity of the rear disk and Vq^ 


is the 


free-stream velocity of the front disk. The power AP added to the 
stream hy any single disk is given hy the change in energy from its 
free-stream to its final-stream conditions. Considering only an 
isentropic flow process this power is given hy 



( 2 ) 


where the subscript n refers to . the nth disk. The total power P for 
all disks is given hy 


P 

m 



(3) 


where t equals total number of disks. It can he shown that 


Z = iAj 

m 2 \ 


3r 



(4) 
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The total power disk-loading coefficient based on the free -stream condi- 
tions of the front disk and the area of the front disk is given by 


P 


c 



which reduces to 


'M h: 


R 


The over -all efficiency is given by 


n = 


-( Y 3k - T °fN 


(5) 


which reduces to 


n = 



(6) 


For isentropic flow the final static pressure and density are the 
same as the initial free-stream values; consequently, equation (5) 
reduces to 



and equation (6) reduces to 


n = 


2 


1 + 


% 

% 


(7) 


( 8 ) 


Equations (7) and (8) permit the calculation of the maximum perform- 
ance for disks in tandem with isentropic flow throughout the stream. For 
these conditions the maximum value of Mg can be 1.0 for an infinite 
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number of disks and the maxlrmj.ni value of 



is the value when the 


flow into the front disk is choked. Under these conditions the maximum 
values of P c and t) then become 


P 


c max 



( 9 ) 


and 


^max 



( 10 ) 
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Figure 1.- Slipstream, stations used in the analysis of single actuator 
diet. Free-stream conditions M* V, p, p, H, and E are given 
at station 0. 
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(a) P c =0.09,- P 0 /P= 0; Pc =0.092. 
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(b)P c =O.I4jP 0 /P=0 j Rc=0.092. 
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(c) P C =0.I4 P 0 /P=0.264 ; P c *=0.092 
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Figure 2.- Graphical solution to obtain thrust balance. Mq = 0.85. 
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Power disk-loading coefficient, P c 


Ob) Mach, number. 
Figure 3*- Continued. 
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Power disk-loading coefficient, P c 

(c) Velocity ratio. 

Figure 3*" Continued. 
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Power disk-loading coefficient, P c 

(&) Density ratio. 


Figure 3.- Continued. 
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Power disk-loading coefficient, P c 


(e) Area ratio. 
Figure 3*- Concluded. 
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Power disk- loading coefficient, P c 

(a) Pressure ratio. 

Figure 4.- Variation of slipstream, flow factors with power loading. 

Mq = 0.70. 
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Power disk-loading coefficient, P c 


♦ 


(b) Mach, number. 
Figure It.- Continued. 
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(c) Velocity ratio. 
Figure 4 .- Continued. 
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Power disk-loading coefficient, P c 


(d) Density ratio. 
Figure 4.- Continued . 
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(e) Area ratio. 
Figure h.~ Concluded. 
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(a) Pressure ratio . 

Figure 5.- Variation of slipstream, flow factors with power loading. 

M q = 0.80. 
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Power disk-loocfing coefficient, Pc 
(c) Velocity ratio. 


Figure 5.- Continued. 
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Figure 5.- Continued. 
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Power disk-loading coefficient, P c 

(e) Area ratio. 

Figure 5.- Concluded. 
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Power disk-loading coefficient,^ 


(a) Pressure ratio. 

Figure 6.- Variation of slipstream flow factors with power loading. 

Mq = 0.85. 
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Power disk- loading coefficient ,P C 
(Id ) Mach, number. 


Figure 6.- Continued. 
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(c) Velocity ratio. 
Figure 6.- Continued. 
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Power disk- loading coefficient, Pc 
(d) Density ratio. 

Figure 6 — Continued. 
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Power disk-loading coefficient, P c 


(a) Pressure ratio. 

Figure 7-~ Variation of slipstream flow factors with. power loading. 

Mq = 0.90. 
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Power disk- loading coefficient ,P C 
(b ) Mach number . 

Figure 7*- Continued. 
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Power disk-loading coefficient^ 

(c) Velocity ratio. 

Figure 7 .- Continued. 
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Power disk-loading coefficient^ 

(d) Density ratio. 
Figure 7.- Continued. 
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Power disk-loading coefficient ,P C 


(a) Pressure ratio. 

Figure 8.- Variation of slipstream flow factors with, power loading. 

Mq = 0.95. 
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Power disk-loading coefficient, P c 

(c) Velocity ratio. 
Figure 8.- Continued. 
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(d) Density ratio. 
Figure 8.- Continued. 
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(e) Area ratio. 
Figure 8.- Concluded. 
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Mach number behind normal shock, M 2n 
Mach number for subsonic -wake solution , Mg 


Figure 9 •- Mach number relation for normal shock and correlation with 
solutions from actuator-disk theory. 
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Figure 10.- Disk, pressure-rise coefficient for incompressible and unchocked 

compressible flows. 








Velocity ratio behind disk, \^/V 0 



Figure 11.- Effect of free-stream Mach, number on velocity ratio behind. 

the disk. Fq/F => 0. 
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Power disk-loading coefficient, P c 


(a) Mq » 0.60, 0.70, and. 0.80. 

Figure 12 .- Comparison of average- velocity ratio through disk far incompressible 
and compressible flows up to choked condition. 
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(b) Mq = 0.80, 0.85, 0.90, and 0.95. 
Figure 12.- Concluded. 
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Figure 13.- Variation of maximum, efficiency with power loading for several 

forward Mach numbers. 
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(a) Pressure ratio. 

Figure 14. - Effect of profile losses on tlae slipstream flow factors for 

subsonic wake. Mq = 0.85. 
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Power disk- loading coefficient, P c 

(b) Mach number. 

Figure 14.- Continued. 
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(c) Velocity ratio. 
Figure 14.- Continued 
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(d) Density ratio. 
Figure 14.- Continued. 
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(e) Area ratio. 
Figure l4.~ Concluded. 


ON 


NA.CA EM L53A.07 



0 .04 08 .12 .16 .20 .24 .28 .32 .36 .40 .44 .48 .52 .56 

Power disk-loading coefficient, P c 


Figure 15.- Effect of profile losses on efficiency. Subsonic veke. 
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